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By the first- principles electronic structure calculations, we find that energetically the most favor- 
able cleaved AFe2As2(001) surface (A=Ba, Sr, or Ca) is A-terminated with a (y/2 x y/2)R45° or 
(1x2) order. The (1x2) ordered structure yields a (1 x 2) dimerized STM image, in agreement 
with the experimental observation. The A atoms are found to diffuse on the surface with a small 
energy barrier so that the cleaving process may destroy the A atoms ordering. At the very low tem- 
peratures this may result in an As-terminated surface with the A atoms in randomly assembling. 
The As-terminated BaFe2As2 surface in orthorhombic phase is (y/2 x y/2)R45° buckled, giving rise 
to a switchable (y/2 x \^2)R45° STM pattern upon varying the applied bias. No any reconstruction 
is found for the other As-terminated surfaces. There are surface states crossing or nearby the Fermi 
energy in the As-terminated and (1x2) A-terminated surfaces. A unified physical picture is thus 
established to help understand the cleaved AFe2As2(001) surfaces. 

PACS numbers: 68.35.B-, 73.20.-r, 74.70.Dd, 68.43.Bc 



Recently great interest has been devoted to the in- 
vestigation of superconduction properties of layered iron 
pnictides upon doping or high pressure pj. Both the neu- 
tron scattering and theoretical studies revealed that the 
undoped parent compounds at atmospheric pressure are 
in either collinear or bi-collinear antiferromagnetic or- 
der below either tetragonal-orthorhombic or tetragonal- 
triclinic structural transition temperature [I, 0, Q, HI- 
The surfaces of high-quality single crystal Ba(Sr)Fe2As2 
have also been probed by Scanning Tunneling Microscopy 
(STM) @, 0, H[. However, the contradictory observa- 
tions on the nature of terminated surfaces have been 
found. On a single crystal of BaFe2As2 cleaved at 20K, 
the STM images observed at 4.3K showed a dominant 
(y/2 x ^2)^45° (denoted by y/2 thereafter) order, which 
was interpreted as due to the ordering of As atoms on 
the terminated surface [6]. The STM measurement on 
another BaFe2As2 single crystal cleaved at 120K also re- 
vealed a y/2 order, but it was interpreted as due to the 
ordering of Ba atoms on the terminated surface 0. In 
contrast, the STM measurement observes both y/2 and 
(1 x 2) orders on a SrFe2As2 single crystal[8]. These re- 
flect there will be rich surface properties distinct from the 
bulk. A thorough investigation on the underlying atomic 
structures and electronic structures of the cleaved sur- 
faces is important to resolve these contradictions. This 
is also important to the understanding of experimental 
data measured by Angle-Resolved Photoemission Spec- 
troscopy (ARPES) which is a surface-sensitive probe. 

In this Letter, we present the results of first-principles 
electronic structure calculations on the (001) surfaces of 
AFe 2 As 2 (denoted as A122 thereafter) (A=Ba, Sr, or Ca). 
We find that the A-terminated surface with either y/2 or 
(1x2) order is energetically most favorable in these mate- 
rials. Cleaving at very low temperature may lift A atoms 



free in a large area, resulting in an ordered As layer to 
expose on the surface. Correspondingly, there are the 
surface electronic states induced crossing or nearby the 
Fermi energy. These surfaces lead to different STM pat- 
terns and affect strongly low energy electronic structures 
of the compounds. 

In our calculations the plane wave basis method was 
used We adopted the generalized gradient 

approximation of Perdew-Burke-Ernzerhof [ll| for the 
exchange-correlation potentials. The ultrasoft pseudopo- 
tentials [12[ were used to model the electron-ion interac- 
tions. After the full convergence test, the kinetic energy 
cut-off and the charge density cut-off of the plane wave 
basis were chosen to be 408 eV and 4082 eV, respectively. 
The Gaussian broadening technique was used and a mesh 
of 16 x 16 k-points were sampled for the (lxl) irreducible 
surface Brillouin-zone integration. We modeled the sur- 
face using a periodically repeated slab of six FeAs layers 
and five or seven A layers plus a vacuum layer of 20 to 
25 A with inversion symmetry through the center of the 
slab. We froze the two innermost middle FeAs layers with 
one A layer and allowed all other atoms to relax. 

In this Letter we call a surface unit cell in reference 
to a square lattice with ~4A periodicity since both As 
atoms and A atoms sit at such a square lattice in the 
bulk. A (lxl) unit cell is thus a 4Ax4A square. After 
cleaving, the surface of AFe2As2 can terminate either at 
an A- atom layer or at an As- atom layer. An Fe-atom 
layer will not expose on the surface due to the strong 
chemical bonding between Fe and As ions. On average, 
half of A atoms should remain on the surface to balance 
the chemical valence. These A atoms can only form ei- 
ther y/2 or (1x2) order, derived by half of A atoms being 
removed alternatively from the original (lxl) square lat- 
tice of A atoms along the square-side or square-diagonal 
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TABLE I: AFe 2 As 2 (001) surface energy differences (A=Ba, 
Sr, or Ca) between the combinations of the two (y/2 x 
V / 2)i?45° A-terminated surfaces (denoted as 2x(y / 2)), the two 
(1x2) A-terminated surfaces (denoted as 2x(lx2)), and the 
(lxl) A- terminated and (lxl) As-terminated surfaces (de- 
noted 2x(lxl)) in nonmagnetic tetragonal phase (tetragonal) 
and antiferromagnetic orthorhombic phase (orthorhombic), 
respectively. Here we set the surface energy of (2x(lxl)) to 
zero. The unit is meV/(l x 1 cell), whereas lxl cell is a 
square. 



AFe 2 As 2 
(001) 


tetragonal 


orthorhombic 


2x(^2) 


2x(lx2) 


2x(lxl) 


2x(>/2) 


2x(lx2) 


2x(lxl) 


BaFe 2 As 2 


-254.8 


-151.0 


0.0 


-325.0 


-221.8 


0.0 


SrFe 2 As 2 


-181.9 


-97.4 


0.0 


-261.9 


-194.7 


0.0 


CaFe 2 As 2 


-101.0 


-107.8 


0.0 


-224.8 


-210.1 


0.0 



direction. In reality, there always yield two surfaces after 
cleaving. It is expected that the cleaved surface would 
be A-terminated with a \pi or (1 x 2) order. On the 
other hand, if one of cleaved surfaces is As-terminated, 
the other will be A-terminated with all the A atoms re- 
maining at (1 x 1) lattice. 

In order to clarify the cleaved surface structure, we cal- 
culated the A-terminated surfaces with (1 x 1), V% and 
(1x2) orders described above and the As-terminated 
surface with (lxl) order, respectively. Table [J sum- 
marizes the numerical results. We find that energeti- 
cally the most favorable cleaved A122 (001) surface is A- 
terminated. For Bal22 it is V^-ordered while for Cal22 
(Srl22) it is mainly (1 x 2) (^-ordered with partly y/2 
((1x2)) order coexisting as the energy difference between 
the (1x2) and y/2 orders becomes small. This may be 
due to the electropositivity becoming smaller in the line 
of Ba, Sr, and Ca. 

Realistically, the kinetics of cleaving process plays an 
important role in the formation of the cleaved surfaces 
of a sample. To understand this more clearly, we carried 
out a calculation to simulate a cleaving process. We di- 
vided a sample into two parts by separating two neighbor 
FeAs layers adiabatically away from each other. The A 
layer between the two neighbor FeAs layers is either di- 
vided half by half into each of the FeAs layers or kept as 
a whole in one of FeAs layers in the separation process. 
It turns out that no significant separation barrier ener- 
gies are found in these two cases. This suggests that the 
energetics shown in Table U may well describe the cleaved 
surfaces. 

We further studied the diffusion of an A atom on the 
As-terminated surface. The energy barrier to move an A 
atom from one equilibrium position to another is found 
to be small so that the fast cleaving process may destroy 
the A atoms ordering. In particular, a very low temper- 
ature cleaving can yield a metastable As-terminated sur- 
face with randomly assembled A atoms. The low energy 
electronic structure is expected to be strongly modified 
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FIG. 1: (Color online) The calculated surface electronic band 
structure of AFe 2 As 2 (A=Ba, Sr, or Ca) (001) surface re- 
ported along the high-symmetry lines of the irreducible sur- 
face Brillouin zone. Fermi level sets to zero. Shaded areas 
correspond to the surface-projected bulk states while thicker 
blue lines correspond to surface states. The surface Brillouin 
zones are given next to the surface band structures whereas 
the dotted line squares are (lxl) surface Brillouin zone, 
(a) (lxl) As-terminated BaFe 2 As 2 surface in the tetragonal 
phase; (b) (y/2 x y/2)R45° As-terminated BaFe 2 As 2 surface 
in the orthorhombic phase; (c) (V2 x y/2)R45° Ba-terminated 
BaFe 2 As 2 surface in the orthorhombic phase; (d) (1x2) Ca- 
terminated CaFe 2 As 2 surface in the tetragonal phase. Note 
that (lxl) surface Brillouin zone is folded in (b), (c), and 

(d). 











FIG. 2: (Color online) Simulated STM images of the (y/2 x 
V / 2)i?45° -ordered As-terminated BaFe2As2 (001) surface in 
orthorhombic phase with bias voltages of (a) 30meV, (b) - 
30meV, and (c) -lOOmeV. The large squares indicate (y/2 x 
y/2)R4b° surface unit cells, while the small square indicates 
(lxl) surface unit cell. Note that the simulated STM im- 
ages for As- terminated Sr(Ca)Fe2As2 surface show the (lxl) 
order, similar to (b). 

on the As-terminated surface since the unbalanced chem- 
ical valence in the surface FeAs layer will induce surface 
states crossing or nearby the Fermi energy. This may 
further induce surface reconstruction through buckling 
or dimerization of the surface As atoms. 

In the tetragonal phase, we don't find any reconstruc- 
tion like buckling or dimerization on the As-terminated 
A122 surfaces except small inwards relaxation. Neverthe- 
less, a number of surface states are induced within the 
bulk energy gap regions, especially those partially filled 
surface states nearby the Fermi energy (see Fig. DJa)). 

In the orthorhombic phase, the As-terminated A122 
surface takes a \/2 magnetic surface unit cell. Our cal- 
culations show that there is neither buckling nor dimer- 
ization except in Bal22 where a small buckling of ~0.lA 
occurs with an energy gain of about 3.9 meV/ (lxl cell). 
Such a small buckling can nevertheless affect strongly 
the STM measurement, as shown in Fig. O For a pos- 
itive bias or a small negative bias (Fig. [2^i) , the calcu- 
lated STM image shows a y/2 order, in which the buck- 
led inward As atoms show strong bright spots while the 
buckled outward As atoms are (nearly) invisible. For a 
mediate bias between -30 meV and -50 meV (Fig. (2Jd), 
a (1 x 1) order of bright buckled As is found, in which 
all the buckled As atoms show the bright spots. For a 
large negative bias below -50 meV (Fig. [2k), a y/2 or- 
der, in which the buckled inward As atoms are nearly 
invisible while the buckled outward As atoms are bright- 
est, is found. Here a positive (negative) bias means the 
STM probes the states above (below) the Fermi energy. 
These results show that the STM image can be switched 
from a \[2 pattern shifted by ~4A to another y/2 pat- 
tern through a (1 x 1) order by varying the bias voltage. 
For the As-terminated surface in Sr(Ca)Fe2As2, the cal- 
culated STM images show a (1 x 1) order, similar to the 
one shown in Fig. [2(b). 

The surface band structure (Fig. H(b)) shows that 
there are a number of surface states induced within the 
bulk energy gap regions in the As-terminated Bal22. 
Here we remind that the nonmagnetic tetragonal surface 




FIG. 3: (Color online) Simulated STM images of the (y/2 x 
y/2)R4b° ordered Ba-terminated BaFe2As2 (001) surface in 
orthorhombic phase with bias voltages of (a) lOOmeV and 
(b) -lOOmeV. The squares indicate (y/2 x y/2)R4b° surface 
unit cells. Note that the similar simulated STM images are 
obtained for the (y/2 x y/2)R45° ordered Sr(Ca)-terminated 
Sr(Ca)Fe2As2 surface. 



Brillouin zone is reduced by half through folding in the 
antiferromagnetic orthorhombic surface Brillouin zone. 
As shown in Fig. [1] (b), the tetragonal M and T become 
T f in the folded Brillouin zone, and the tetragonal X is 
now located at M' . 

The buckling of a surface is to gain energy by lifting the 
degeneracy of energy band due to breaking of mirror sym- 
metry. The buckling of pairs of As atoms at the surface 
causes the surface atoms to partially dehybridize. This 
causes some electronic charge transfer from the inward 
atoms to the outward atoms so that the initial partially 
filled surface states become insulating (see Fig. [1(b)). 
This is the reason why no buckling happens in the non- 
magnetic tetragonal phase whereas the energy gain by 
the charge transfer is smaller due to the stronger metal- 
licity. This is also the reason why the STM can observes 
the inward and outward atoms at a positive and negative 
bias, respectively. On the other hand, the buckling will 
cause a loss of the surface lattice distortion energy. These 
competing effect determines whether or not the buckling 
will take place. In contrast, for Sr(Ca)122, no recon- 
struction is found in the As-terminated surface. This is 
because in Sr(Ca)122 the lattice constants are smaller 
than that in Bal22 and that the lattice distortion will 
cost more energy than the buckling energy gain. This 
is similar to the buckled Si (001) surface versus the non- 
buckled C(001) surface. As a result, the As-terminated 
Sr(Ca)122 surfaces always yield a (1 x 1) order in either 
tetragonal or orthorhombic phase, similar to Fig. [2(b). 

For the y/2- ^-terminated A122 surface, no lattice re- 
construction (buckling or dimerization) is found. One 
expects that those surface states crossing the Fermi en- 
ergy in the As-terminated surfaces (see Fig. [T] (a) and 
(b)), induced by the unbalanced chemical valence, would 
be eliminated by the A atoms donating charge into the 
FeAs layer to make the chemical valence balanced in the 
whole surface layer. This is verified by our calculations, 
as shown in Fig. [T] (c). The remaining surface states are 
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FIG. 4: (Color online) Simulated STM images of the (1 x 2) 
ordered Ca-terminated CaFe2As2 (001) surface in orthorhom- 
bic phase with bias voltages of (a) lOOmeV and (b) -lOOmeV. 
The squares indicate (1x2) surface unit cells. Note that the 
similar simulated STM images are obtained for the (1x2) 
ordered Sr(Ba)-terminated Sr(Ba)Fe2As2 surface. 

mainly distributed around the edges of the bulk energy 
gap regions, induced by the reduced surface symmetry. 
The calculated STM images show a V^-c-rdering pattern 
in both tetragonal and orthorhombic phases (see Fig|3|), 
in agreement with the recent STM observation [7j. 

For the (1 x 2)-A-terminated A122 surfaces, no any 
other reconstruction (buckling or dimerization) is found 
for the surface A atoms either. However, in both tetrag- 
onal and orthorhombic phases, the surface A atoms will 
move downwards so that the underneath As atoms move 
closely with each other by ~0.05 A. This small dimer- 
ization is perpendicular to the surface A atom row and 
located in the troughs between the A atom rows, and lead 
to some new surface states near the Fermi energy (Fig. 
DJd)). However, the calculated STM images show the ap- 
parent dimer rows along the surface A atom rows rather 
than the troughs, due to the hybridization between the 
surface A atoms and As atoms (Fig. [4]), also in agreement 
with the recent experimental STM observation [8j. 

In summary, a unified physical picture on cleaved 
AFe2As2 (001) surfaces (A=Ba, Sr, or Ca) is estab- 
lished based on the first principles band structure calcu- 
lations. The energetically most favorable cleaved surface 
is A-terminated with (y/2 x y/2)R45° or (1 x 2) order. 
For Bal22, the surface A atoms are found to be in the 
{y/2 x \^2)R4:5 order, in agreement with the experimen- 
tal observation [7]. However, as there is a small diffuse 
energy barrier for A atoms, a fast cleaving at very low 
temperature may yield an As-terminated ordered surface 



with randomly assembled A atoms as a metastable sur- 
face structure. This As-terminated surface is also found 
to be (V2xV2)R45° ordered. This leads to a natural ex- 
planation to the surface structure of Bal22 cleaved at low 
temperature measured by Pan and coworkers [6j . Further- 
more, we predict a switchable (y/2 x ^2)^45° STM pat- 
tern upon varying an applied bias, which can be a hall- 
mark for testing whether the As-terminated or the Ba- 
terminated is. This theoretical prediction can be further 
verified by annealing an As-terminated surface to high 
temperatures so that the A atoms can diffuse to form an 
^-terminated surface. For the As-terminated Sr(Ca)122 
surface, there is no buckling, dimerization, or any other 
reconstruction. Thus the STM-observed coexistence of 
(y/2 x V2)R45° and (1 x 2) patterns on Srl22 surface^ 
should be attributed to the Sr-termination. For Cal22, 
we find that the Ca-terminated surface is (1 x 2) ordered 
which coexists partly with a (y/2 x \/ r 2)R4b° order. Fi- 
nally, among these surface structures, there are essential 
surface states induced crossing or nearby the Fermi en- 
ergy within the bulk energy band gap regions, which will 
affect those surface sensitive experimental probing. 
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